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A  Measure  of  Rotatability  for 


Response  Surface  Designs 


A.  1.  Khuri 


Department  of  Statistics 
University  of  Florida 
Gainesville,  FL  32611 


ABSTRACT 

A  measure  is  introduced  in  this  paper  which  quantifies  the 
amount  of  rotatability  in  a  given  response  surface  design.  The 
measure,  which  is  expressible  as  a  percentage,  takes  the  value  100 
if  and  only  if  the  design  is  rotatable.  One  of  the  main 
advantages  of  this  measure  is  that  it  can  be  used  to  ''repair"  a 
nonrotatable  design  by  the  addition  of  experimental  runs  which 
maximize  the  percent  rotatability  over  a  spherical  region  of 
interest.  Four  numerical  examples  are  given  to  illustrate  the 
applications  of  this  measure. 


KEY  WORDS:  Nonrotatable  designs;  Percent  rotatability;  Repairing 
rotatability;  Design  moments;  Design  augmentation; 
Cone  of  rotatability. 
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1.  INTRODUCTION 


Consider  fitting  a  linear  response  surface  model  of  order  d 
in  k  input  variables,  ,X2 , . • • ,x^,  over  a  spherical  region  of 
interest,  R,  using  a  design  consisting  of  n  experimental  runs. 

This  model  can  be  written  in  vector  form  as 

E(^)  =  XJ3.  0-0 

where  ^  is  a  vector  of  n  observations,  E(^)  denotes  the  mean  or 
expected  value  of  X  is  an  nxp  matrix  of  rank  p  whose  elements 
are  known  functions  of  the  design  settings  of  the  input  variables, 
and  J3  is  a  vector  of  unknown  regression  coefficients.  We  assume 
that  the  variance-covariance  matrix  of  j  is  given  by  o2  1^ ,  where 
a2  is  unknown  and  I  is  the  identity  matrix  of  order  n*n.  We 
denote  by  D  the  nxk  design  matrix  whose  u1'1  row  consists  of  the 
settings  of  the  k  input  variables  at  the  u^  experimental  run 
(  u— 1,2 . n) . 

The  predicted  response  value  at  a  particular  point  x  = 

(xj ,X2 , . • • jX^) '  in  the  region  R  will  be  denoted  by  y(x).  This 
value  is  obtained  by  substituting  the  elements  of  8  in  model  (1.1) 
by  the  corresponding  elements  of  8,  the  least  squares  estimator 
of  8,  namely, 

8  =  •  (1.2) 
The  design  J3  is  said  to  be  rotatable  if  the  variance  of 
y(x)  ,  which  in  general  depends  on  x , ,Xj , • • • ,x,  ,  is  a  function  of 


only  the  distance  of  the  point  £  from  the  center  of  the  design. 
Thus,  when  the  design  is  rotatable,  the  prediction  variance, 
denoted  by  var[y(x)],  is  the  same  at  all  points  x  which  are  equi¬ 
distant  from  the  design  center.  Consequently,  in  the  space  of  the 
input  variables,  surfaces  of  constant  prediction  variance  form 
concentric  hyperspheres  (circles  and  spheres  in  two-dimensional 
and  three-dimensional  Euclidean  spaces,  respectively). 

The  concept  of  rotatability  was  first  introduced  by  Box  and 
Hunter  (1957)  and  has  since  become  an  important  design  criterion. 
One  of  the  desirable  features  of  rotatability  is  that  the  quality 

A 

of  prediction,  as  measured  by  the  size  of  var[y(x)],  is  invariant 
to  any  rotation  of  the  coordinate  axes  in  the  space  of  the  input 
variables.  Furthermore,  if  optimization  of  y(a)  is  desired  over 
concentric  hyperspheres  within  a  certain  region  of  interest,  it 
would  be  very  desirable  to  have  a  rotatable  design.  Otherwise, 
poor  estimates  of  the  optimum  might  result  (see  Khuri  and  Myers, 
1979) 

A  simple  characterization  of  rotatability  is  given  in  terms 
of  the  elements  of  the  X'X  matrix.  We  shall  refer  to  these  ele- 

<V  <~v» 

raents  as  design  moments  (traditionally,  the  elements  of  the 
matrix  X'X/n  are  referred  to  as  design  moments).  In  general,  a 
design  moment  for  a  model  such  as  (1.1)  of  order  d  and  in  k  input 

SiSt  <5 
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where  6^  •  •  •  ,6^  are  nonnegative  integers  and  xuj  is  the  level 

of  the  jC^  input  variable  used  in  the  uC^  experimental  run 

k 

( j  =  l , 2 . k;  u=  1 , 2  , . . . ,n) .  The  sum,  Z  6.,  is  called  the  order 

j  =  l  J 

of  the  design  moment  and  is  denoted  by  6(6  =  0,1,..., 2d).  For 

,  2  3 

example,  ( 1  3  5  )  is  a  design  moment  of  order  6=6  and  is  equal 


to  Z  x  ,x  ,x  .. 
.  ul  uJ  u5 
u=l 


A  necessary  and  sufficient  condition  for  a  design  for  fitting 
a  model  such  as  (1.1)  (of  order  d  and  in  k  input  variables)  to  be 
rotatable  is  that  the  design  moments  of  order  6 (6=0 , 1 , .  . .  , 2d )  be 
of  the  form 


(1 


6i  6 
2 


0, 


0 r  .n. 

-A  Jr.1.. 


2S/2  .n 
.1  =  1 


6  j  ! 


tV2)! 


if  any  6^  is  odd 

(1.4) 

if  all  of  the  6 . 's  are  even, 

t  9 


where  6^  is  a  quantity  which  depends  on  d,6,  and  n  (see  Box  and 
Hunter  1957,  Myers  1976,  ch.  7).  For  convenience,  we  say  that  a 
design  moment  is  odd  if  at  least  one  of  its  6j's  is  odd  and  a 
design  moment  is  even  if  all  of  the  6 j ’ s  are  even.  Note  that  a 
design  moment  of  order  6  =  0  is  equal  to  n.  A  design  whose 


moments  do  not  conform  to  formula  (1.4)  is  said  to  be  nonrotat- 


Quite  often,  a  nonrotatable  design  may  exhibit  surfaces  of 
constant  prediction  variance  which  are  nearly  spherical.  In  this 
case,  the  design  is  described  s  being  near  rotatable.  This 
occurs,  for  example,  when  a  rotatable  design  is  deformed  due  to 
incorrect  settings  of  some  of  the  input  variables,  or  because 
certain  specified  levels  of  the  input  variables  may  be  difficult 
to  employ  in  practice.  In  another  situation,  a  rotatable  design 
may  undergo  certain  modifications  to  fit  the  needs  of  the  experi¬ 
ment.  The  modifications  might  involve  adding  new  design  points, 
or  shifting  existing  design  points,  in  order  to  gain  more  informa¬ 
tion  in  a  certain  region  of  interest  (see  Littell  and  Mott  1974). 
In  such  situations  it  is  of  interest  to  assess  the  effects  of 
deformation  or  modification  on  the  property  of  rotatability. 

To  assess  the  degree  of  rotatability,  it  has  been  customary 
to  inspect  contour  plots  (in  case  of  k=2  input  variables)  of 
constant  prediction  variance  to  see  how  close  they  are  to  being 
circular.  Such  a  practice  .however,  has  its  limitations  when  the 
number  of  input  variables  exceeds  3  in  addition  to  being  subjec¬ 
tive.  In  this  paper,  we  provide  a  quantitative  measure  of  rotat¬ 
ability  for  a  response  surface  design.  This  measure  takes  values 
between  0  and  100  with  the  latter  value  being  attained  when  the 
design  is  rotatable.  The  proposed  measure  can  be  useful  in  the 
following  situations: 

(i)  to  quantify  the  degree  of  rotatability  of  a  nonrotatable 


A1! 
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design  so  that  a  determination  of  how  close  the  design  is 
to  being  rotatable  can  be  made. 

(ii)  to  compare  designs  on  the  basis  of  their  degrees  of 
rotatability. 

(iii)  to  assess  the  extent  of  departure  from  rotatability  when  an 
already  rotatable  design  is  deformed  or  modified. 

(iv)  to  repair  a  nonrotatable  design. 


2.  A  MEASURE  OF  ROTATABILITY 

Consider  again  the  model  given  in  (1.1),  which  if  we  recall 
is  of  order  d  in  the  input  variables  Xj  ,X2>  • « •  »xjc.  Let  us  suppose 
that  in  the  spherical  region  R  over  which  this  model  is  fitted  the 
input  variables  have  been  coded  so  that 


E  z  .  =  0, 

i  UJ 
u=l  J 


j-l,2,...,k 


j-l,2,.«. , k, 


(2.1) 


where  z^j  denotes  the  coded  value  of  xyj ,  the  actual  value  of  the 

uth  level  of  variable  j,  and  a  is  some  positive  constant.  The 

coding  can  be  accomplished  by  applying  the  transformation,  z  .  = 

n 

(x  x. )/s • ,  where  x.  =  lx  . /n  and  s.  is  given  by 

J  J  J  J  i  la  I  J 

u=l  J 


II 

ij  =  [  E  (x  x  . 
J  u=l  UJ  J 


By  adapting  the  notation  described  in  (1.3)  to  the  coded  variables 


the  equalities  described  in  (2.1)  can  be  written  as 


(j)  =  0  , 

j  —  l,2,...,k. 

(j2)  =  a  , 

j  =  1,2,.. . , k. 

Under  this  coding  scheme,  the  center  of  the  design  coincides  with 
the  point  at  the  origin  of  the  coordinates  system.  The  coding 
also  helps  to  standardize  the  input  variables  which  may  have 
different  units  of  measurements..  Furthermore,  the  spread  of  the 
design  will  be  the  same  in  all  directions  of  the  coordinate  axes. 

In  terms  of  the  coded  variables,  model  (1.1)  can  be  expressed 
as 

E(*>  =  a-  <2.3) 

where  the  matrix  jg  is  of  the  same  form  as  X  except  that  the  zuj's 
are  used  instead  of  the  x^'s,  and  y  is  the  new  regression  coeffi¬ 
cient  vector.  Throughout  the  remainder  of  this  paper,  the  design 
moments  defined  in  (1.3)  will  be  formulated  in  terras  of  the  coded 
variables . 

We  shall  now  introduce  a  measure  of  rotatability  for  a  given 
response  surface  design.  Such  a  measure  should  be 

(i)  a  function  of  the  levels  of  the  input  variables  used  in  the 
design, 

(ii)  invariant  to  the  value  of  the  scale  parameter,  a,  in  formula 

(2.1) 

(iii)  invariant  to  the  addition  of  points  at  the  center  of  the 
design  (which  coincides  with  the  origin  of  the  coordinates 


system  by  the  coding  scheme). 

The  reason  for  condition  (ii)  is  that  there  can  be  an  infinite 
number  of  designs  that  are  derivative  of  one  another  by  changing 
the  value  of  a.  Since  such  designs  are  not  different  with  respect 
to  rotatability,  their  measures  of  rotatability  should  be  the 
same.  As  for  condition  (iii),  rotatability  depends  only  on  design 
moments  of  order  6  >  1  and  these  remain  invariant  to  the  addition 
of  center  points. 

Three  major  steps  are  involved  in  the  development  of  the 
proposed  measure  of  rotatability.  Details  of  these  three  steps 
are  given  below. 

Step  1:  The  Reduction  of  All  Design  Moments  to  Scale-Free 

Quantities 

Let  us  denote  by  v(Z'Z)  the  vector  consisting  of  the  elements 
of  Z'Z  which  are  located  above  the  diagonal  as  well  as  those  along 
the  diagonal.  This  vector  is,  therefore,  of  dimension 
P*  =  p(p+l)/2,  where  p  is  the  number  of  parameters  in  model 
(2.3).  Its  elements  are  obtained  by  listing  the  elements  in  the 
first  row  followed  by  those  in  the  second  row,  etc.,  of  the  upper 
triangular  half  of  Z'Z.  It  can  be  shown  that  the  element  at  the 
(i,j)t'n  location  (j  >  i)  in  this  upper  triangular  half  is  the  i. 
element  of  the  vector  v(Z'Z),  where  i  =  f(i,j)  is  the  function 


S.  S 


f(i,j)  =  (i-1 ) [p-( i/2 )  ]  +  j,  j  >  i.  (2.4) 

For  example,  for  a  second-order  model  of  the  form 

E(y)  -  V  Vl+  y2z2*  y[2v2+  Tr„^+  y22z\, 
p  =  6,  v(Z'Z)  is  of  order  21x1,  and  the  (3,5)*"^  element  of  Z^Z, 
which  is  equal  to  the  design  moment  (1  2),  is  the  14^  element 
of  v(Z'Z). 

<v  /w  »v 

We  note  that  the  vector  v(Z'Z)  is  scale  dependent  since  it 

(v  «v 

depends  on  the  value  of  a  used  in  formula  (2.1).  To  remove  this 
dependency,  we  divide  each  design  moment  of  order  6(6=  0,1,..., 2d) 


by  t  ,  where  t  is  given  by 


t  =  [  Z  (j2)/k] 

j  =  l 


1/2 


(2.3) 


that  is,  t2  is  the  average  of  all  even  design  moments  of  order 
2.  We  note  that  by  the  coding  scheme  described  in  (2.2),  i  = 
a  '  .  This  operation  amounts  to  premultiplying  v(Z'Z)  by  a 
diagonal  matrix,  denoted  by  £,  °f  order  p*  x  p*.  The  diagonal 
elements  of  A  are  in  a  one-to-one  correspondence  with  the  dements 
of  v(Z'Z):  If  an  element  of  v(Z'Z)  is  a  design  moment  of  order 
6(6=0,  1, ... ,2d) ,  then  the  corresponding  diagonal  element  o!  £  is 
equal  to  1/t  .  Hence,  the  elements  of  the  vector 

u(Z'Z)  =  A  v(Z’Z)  (2.6) 

are  scale-free  quantities. 
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Step  2:  The  Introduction  of  a  Canonical  Representation  of  the 

Z'Z  Matrix  for  a  Rotatable  Design 

If  the  design  for  model  (2.3)  is  rotatable,  then  its  design 
moments  of  order  6(6=  0,1,..., 2d)  must  have  the  special  form 
described  in  (1.4).  Hence,  all  the  elements  of  v(Z'Z)  corre- 

<V  <v  /V 

sponding  to  odd  design  moments  should  be  equal  to  zero  and  all 
those  corresponding  to  even  design  moments  must  have  the  form 
0^c(6j  ,6  , ... ,6^),  where 

k 

.  n  6  .  !  k 

c(6  ,6  ,...,6  )  =  -1 1  •--J -  ,  E  6 .  =  6  (6:-0, 2  , . . .  ,  2d ) .  (2.7 

2 6 / 2  J  (6./?)!  j=l  J 
j=  1  J 

Let  Zr  denote  the  7j  matrix  when  the  design  is  rotatable.  In  this 

case  the  vector  v(Z'X  )  can  be  represented  as 
~  ~r~r 


&&  =  °oV  82^2  + 


+  6  2d~  2d  ’ 


(2.8) 


where  w  (6=0 , 2  , . . .  ,  2d)  is  a  vector  of  order  p** 1  whose  elements 

are  in  a  one-to-one  correspondence  with  the  elements  of  v(Z'Z  ): 

~  ~  rv'  r 

All  those  corresponding  to  design  moments  of  order  different  from 
6  and  odd  design  moments  of  order  6  are  equal  to  zero,  whereas 
those  corresponding  to  even  design  moments  of  order  6  have  values 
given  by  formula  (2.7).  In  other  words,  depends  only  on  the 
values  of  the  6 j ’ s  which  designate  the  even  design  moments  of 
order  6 ( 6 =0 , 2 , . . . , 2d )  of  a  rotatable  design.  From  (2.7)  we  note 


that  the  elements  of  w  are  all  zero  except  for  the  first  element 


which  is  equal  to  one.  Also,  from  (1.4)  it  can  be  seen  that  0q 
n  and  that  0£  =  a  by  the  coding  scheme  of  formula  (2.2).  Hence, 
by  (2.5),  @2  =  t2  . 

Since  the  elements  of  the  diagonal  matrix  $  in  (2.6)  that 

correspond  to  design  moments  of  order  6  are  equal  to  l/t1^,  then 

from  (2.6)  and  (2.8)  the  following  canonical  representation 

of  u(Z'Z  )  can  be  obtained: 

~~r~r 


u(Z'Z  ) 
~  ~r~r 


d 

Z  0 


m=0 


2m  ~2m/T 


2m 


=  n  o>  +  u)  +  E  k  u>  , 
~0  ~2  „  2m  ~2m 

m=2 


(2.S 


ry 

where  <2m  =  ®2m^TZm’  or  e9uivalently> 

k 2m  ~  ®2m^®2  »m=2,3,...,d.  (2. 1C 

Since  02m  is  a  design  moment  of  order  5=  2m  as  can  be  seen  from 

(1.4)  by  taking  m  of  the  <5j's  equal  to  2  and  the  remaining  6j's 

equal  to  zero,  the  <2m's  (2.10)  are,  therefore,  scale  free. 

The  parameters  K4 .  •  •  •  »*2d  can  c^osen  by  the  experimenter 

depending  on  whatever  additional  properties  the  rotatable  design 

is  required  to  have.  It  is  to  be  noted  that  the  's 

~2m 

(m=0, 1 , . . . ,d)  in  (2.9)  are  pairwise  orthogonal,  that  is, 

oil  a>„  ,  =  0  for  m  *  m'.  This  follows  from  the  fact  that  all  the 
~2m  ~2m 

elements  of  are  zero  except  for  those  corresponding  to  even 

design  moments  of  order  2m.  For  these  latter  elements,  the 


corresponding  elements  of  0)?  , , 


for  m*m',  must  be  zero  by 


definition.  Consequently,  the  vectors  w  ,  w  w  span  a 

~4  ~2d 

(d-1 )-dimensional  EucLidean  space.  Furthermore,  since  <2m  >  0  for 
in  2 , 3 , . . .  ,d ,  then 


V  =  E  K  W 
~  ~  2m  ~2m 

171  =  2 


(2.11) 


represents  a  vector  in  a  closed  convex  cone  K  in  this  Euclidean 
space.  By  definition,  a  closed  subset  S  of  a  Euclidean  space  is  a 
closed  convex  cone  if  for  any  vectors,  Xj  and  x^,  in  S  and  any 
nonnegative  scalars,  A^  and  A2,  the  vector  A^x^  +  ^  2~2  belongs  to 
S.  We  shall  refer  to  the  cone  K  as  the  cone  of  rotatabi lit y . 


Step  3:  The  Derivation  of  the  Rotatability  Measure 

Let  us  now  suppose  that  we  have  an  n-point  design,  D,  not 
necessarily  rotatable,  for  fitting  model  (2.3).  As  before,  the 
input  variables  are  coded  as  in  (2.2).  The  corresponding 
vector  u(Z'Z)  in  (2.6)  can  then  be  written  as 

u(Z'Z)  =  n  w0+  w2+  u*( Z '  Z) ,  (2.12) 

where  u.  and  are  the  same  as  in  (2.9),  which  holds  for  a  rotat- 
able  design  having  the  same  number,  n,  of  runs.  The  elements 
of  u*(Z'Z)  are  equal  to  the  corresponding  elements  of  u(Z'Z) 
except  for  those  that  correspond  to  design  moments  of  order  6=0 
and  even  design  moments  of  order  6  =  2,  which  are  equal  to  zero. 

To  measure  the  rotatability  of  t  fie  design  JJ  it  is  necessary 
to  determine  how  well  the  vector  u*(Z'Z)  can  he  approximated  with 


a  vector  v  of  the  form  given  in  (2.11).  This  is  equivalent  to 

finding  a  vector  v  e  K,  the  cone  of  rotatability  described 

earlier,  that  is  closest  (in  terms  of  the  Euclidean  norm) 

to  u*(Z'Z).  For  this  purpose  the  parameters  k, ,  k„.  in 

~~~  r  4  6  2d 

(2.11)  are  chosen  so  that 

d 

Q  (D)=  H  u*(  Z '  Z)  -  £  k„  in.  H2  (2.13) 

n  ~  ~  ~  ~  „  2m  ~2m 

m=2 


is  minimum,  where  11*11  denotes  the  Euclidean  norm,  or  length,  of  a 
vector.  In  Appendix  A  it  is  shown  that  the  minimum  value  of  Qn(J?) 
is  given  by  the  formula 


min[Q  (D)]  =  Hu*(Z’Z)ll2  -  £  [u* 1  (Z'Z)w2m]2 /Doj^ll2  .  (2.14) 

ra=2 


It  is  interesting  to  note  that  with  being  given  as  in  (A. 3)  in 

d 

Appendix  A,  the  vector  £  K2m^2m  Pr °-3 e c ^ i °n  of  the 


m=2 


vector  u*(Z'Z)  on  the  cone  of  rotatability  K.  The  square  of  the 
Euclidean  norm  of  this  projection  is  given  by  the  absolute  value 
of  the  second  term  on  the  right  side  of  (2.J4)  and  represents  the 
portion  of  llu*(Z'Z)ll2  which  can  be  attributed  to  rotatability. 
Hence,  as  a  measure  of  rotatability  for  the  design  D  I  choose  the 
quantity 


$n(D)  =  100{llu*(Z’Z)ll2  -  min[Qn(D)]}/iiu*(Z'Z)ll2 
d 

=  100{  £  [u*'(Z'Z)w2m]2/lla;2mn2}/»u*(Z'Z)ll2  ,  (2.15) 

m=2 


which  represents  the  percent  contribution  of  rotatability  to  the 


magnitude  of  llu*(Z'Z)ll2 .  In  other  words,  $n(]D)  represents  the 
percent  rotatability  that  is  inherent  in  the  design  J).  This  is 
analogous  to  the  use  of  the  coefficient  of  determination,  R2 ,  in 
regression,  which  measures  the  proportion  of  the  total  variation 
in  the  response  that  is  caused  by,  or  attributed  to,  the  fitted 
model.  If  jD  is  rotatable,  then  u*(Z'Z)  must  belong  to  the  cone  of 
rotatability  K,  hence  min[Qn(D)]  =  0  and  $n(D)  =  100.  A  large 
value  of  $nQD)  is,  therefore,  an  indication  that  the  design  JD  is 
near  rotatable. 

The  measure  of  rotatability  defined  in  (2.15)  satisfies 
conditions  (i),  (ii),  and  (iii)  stated  earlier  in  this  section. 
This  follows  from  the  fact  that  the  vector  u*(Z'Z)  is  scale  free 

~  rn  /v 

and  its  elements  depend  on  design  moments  of  order  two  or  higher, 
hence  they  are  unaffected  by  the  addition  of  center  points. 

3.  REPAIRING  ROTATABILITY 

The  measure  developed  in  Section  2  can  be  used  effectively  to 
increase  the  percent  rotatability  of  a  nonrotatable  design  by 
augmenting  it  with  additional  runs  chosen  appropriately.  This  is 
particularly  useful  in  situations  where,  because  of  technical 
limitations,  some  of  the  input  variables  of  a  rotatable  design  are 
set  at  levels  different  from  those  in  the  original  design  causing 
it  to  become  nonrotatable.  Also,  it  frequently  happens  that  a 
rotatable  design  is  purposely  modified  by  the  introduction  of  new 


runs  in  order  to  concentrate  information  in  certain  areas  of 


interest  (see  Littell  and  Mott  1974).  This  modification  usually 
results  in  loss  of  rotatability. 

The  choice  of  points  to  be  added  to  a  nonrotatable  design  to 
increase  its  percent  rotatability  can  be  accomplished  in  a 
sequential  manner  as  follows:  Let  denote  a  given 


nonrotatable  design  consisting  of  Oq  experimental  runs.  For  any 
point,  in  the  experimental  region  R,  the  design  consists 
of  the  design  augmented  with  x.  The  percent  rotatability 


of  is  given  by  $  A  new  design  point,  ,  is  now 

Q  ~ 

added  to  to  obtain  the  design  .  This  new  point  is 

(0) 

chosen  by  maximizing  the  percent  rotatability  of  with  respect 

to  x  over  R.  For  simplicity  the  design  is  written  as 

We  thus  have 


$  =  max[$  rD(0))]. 

V1  ~  xeR  n0+1 


A  second  design  point,  #2»  is  subsequently  added  to  '  by 
repeating  the  same  process  as  above,  but  with  replacing 

Let  be  the  design  consisting  of  ^  augmented  with 

£2*  By  continuing  this  process  we  obtain  the  sequence  of  designs, 
^  2^  , . . . , \ . . .  ,  where  is  obtained  from  D^1  ^  by 

augmenting  it  with  x^  (i=l,2,...). 

(  -j  \ 

The  percent  rotatability  associated  with  design  D  is 

$  (DU)),  where  =  nQ  +  i  (i=0,l,...).  In  Appendix  B  it  is 

ni 


shown  that 


4>  (D(l))  <  4>  (D(l  +  1)),  i=0, 1 , . . .  (3.1) 

i  l+l 

Inequality  (3.1)  implies  that  the  sequence,  (Dv  J } ^ _q »  is 

i 

monotonically  increasing.  Since  this  sequence  is  bounded  by  100 

as  can  be  seen  from  (2.15),  it  must  coverge  to  100  (see,  for 

example,  Rudin  1976,  p.  55),  that  is, 

lim  4>  fD(l))  =  100. 
i+°°  i 

In  Khuri  (1985)  it  is  further  shown  that  (3.1)  is  in  fact  a  strict 
inequality. 

The  addition  of  new  runs  to  increase  the  percent  rotatability 
is  tantamount  to  "repairing”  rotatability.  Hobble  and  Mitchell 
(1972)  have  indirectly  done  so  by  maximizing  the  determinant 
| X * X )  through  design  augmentation,  where  X  is  the  matrix  given  in 
(1.1).  The  reason  rotatability  can  be  partially  restored  by  their 
method  is  based  on  theoretical  results  by  Wynn  (1970)  and  Kiefer 
(1961).  More  specifically,  the  method  of  Hebble  and  Mitchell 
leads  in  the  limit  to  a  D-optimal  design  as  was  shown  by  Wynn 
(1970).  On  the  other  hand,  Kiefer  (1961)  showed  that  D-optimal 
designs  over  a  spherical  region  are  rotatable. 

4.  EXAMPLES 

Four  numerical  examples,  4.1  -  4.4,  are  presented  in  this 


section  in  order  to  illustrate  the  applications  of  the  measure  of 


rotatability.  The  purpose  of  example  4.1  is  to  demonstrate  the 
actual  implementation  of  formula  (2.15).  Example  4.2  shows  how 
the  measure  can  be  used  to  compare  designs  on  the  basis  of  their 
degrees  of  rotatability.  Examples  4.3  and  4.4  serve  to  highlight 
the  utility  of  the  measure  in  repairing  rotatabilil  y  with  the 
latter  example  describing  an  actual  experimental  situation. 

Example  4.1 

Consider  fitting  a  second-order  model  of  the  form 

E(y)  -  e0+  BlV  e2x2+  B12x1x2.  6,^2+  822x/  (4.1) 

using  a  32  factorial  design  with  the  jth  factor  having  the  levels 
-1 , 0, 1 ( j= l , 2) .  The  coding  described  in  (2.?)  is  already  satisfied 
with  (j2)  =  6(j=l,2).  In  this  case  a=6,  d=2,  k=2,  and  n=9.  The 
elements  of  the  vectors  v(Z'Z),  u(Z'Z),  u*(Z'Z),  u..  w ,  , 
and  the  diagonal  elements  of  the  diagonal  matrix  A  (see  2.6  and 
2.12)  are  given  in  Table  1.  The  value  of  min[Q  (D)]  in  (2.14)  is 
.00555  and  ^n(D)  in  formula  (2.15)  is  thus  equal  to  93.08.  We 
conclude  that  the  32  factorial  design  is  93.08%  rotatable. 

Example  4,2:  Roquemore's  (1976)  Hybrid  Designs 

These  designs  were  introduced  by  Roquemore  (19/6)  to  emulate 


certain  characteristics  of  central  composite  or  regular  polyhedral 
designs.  They  are  economical  and  are  supposed  to  be  near  rotat¬ 
able.  Roquemore  constructed  three  3-variable  designs,  labeled  by 
him  as  310,  311A,  and  31  IB,  for  a  second-order  model.  These 
designs  are  reproduced  in  Table  2.  He  pointed  out  that  the  31 1A 
design  was  the  most  nearly  rotatable.  This  was  based  on  observing 
contour  plots  of  prediction  variance  and  on  a  comparison  of  the 
values  of  the  ratio  ( j1* )/( j2£2 ) ,  j  *£,  with  the  value  3  (this 
ratio  must  be  equal  to  3  in  order  for  a  design  for  a  second-order 
model  to  be  rotatable  as  can  be  verified  from  1.4).  By  applying 
the  measure  of  rotatability  given  in  (2.15)  it  was  found  out  that 
the  percent  rotatability  values  for  the  310,  31 1A,  and  31  IB 
designs  were,  respectively,  94.89,  99.40,  and  98.99.  This  con¬ 
firms  Roquemore's  assertion  concerning  the  31 1A  design.  It  can 
also  be  seen  that  all  three  designs  have  high  percent  rotatability 
and  that  the  31  IB  design  is  as  nearly  rotatable  as  the  31 1A 
design.  Design  31  IB  was  also  reported  as  the  most  efficient  of 
the  three  designs  with  respect  to  both  the  D-  and  G-optimality 
criteria. 

This  example  clearly  shows  that  in  conjuction  with  other 
measures  of  design  efficiency,  the  rotatability  measure  can  be 
utilized  to  select  a  design  that  possesses  several  characteristics 
of  interest  to  the  experimenter. 
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Example  4.3:  Hebble  and  Mitchell  (1972,  Example  1,  p.  769) 

A  second-order  model  in  two  input  variables  is  fitted  using 
the  ten-point  design  described  in  Table  3.  This  design  was 

originally  planned  as  a  rotatable  central  composite  design  with 
two  center  points.  The  region  of  interest,  R,  is  circular  with 
center  at  (0,0)  and  radius  =  2.  In  terms  of  the  coded  variables, 
Zj  and  Z£,  the  model  is  written  as 

E(y>  -  Y0+  Y,V  V2+  Y12zlV  Y  1 IZ1  +  *22*2  ‘ 


The  percent  rotatability  of  '  according  to  formula  (2.15) 

is  80.65.  This  design  was  subsequently  augmented  with  three 

additional  repair  points,  one  at  a  time,  to  get  the  sequence  of 

designs,  and  D^\  As  was  described  in  Section  3,  the 

design  is  obtained  from  the  design  ^  by  the  addition  of 

the  point  x.  that  maximizes  $  (d^  with  respect  to  £  over  the 

i  ~ 

region  R(i=l,2,3).  The  maximization  of  <*>  was  carried  out  by  using 
SEARCH,  a  computer  program  written  by  Conlon  (1979)  and  is  based 
on  the  controlled  random  search  procedure  introduced  by  Price 
(1977).  The  optimal  locations  of  the  additional  three  points  and 
the  corresponding  percent  rotatability  of  the  ensuing  designs, 

and  JQ^',  are  given  in  Table  4.  Note  that  the  first 


two  points  cause  a  sizable  increase  in  the  percent  rotatability 
beyond  the  value  80.65  for  the  intitial  design  ^  total  of 
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less  than  six  seconds  of  CPU  time  was  needed  to  locate  the 


additional  three  points. 

For  each  of  the  four  designs,  and 

contour  plots  of  var[y( <z)]/o2  were  drawn  within  a  region  in  the 
space  of  the  coded  variables  which  corresponds  to  R,  where  a2  is 
the  random  error  variance  and  z  =  (z^,z^)’.  These  plots  are  shown 
in  Figure  1.  It  can  be  clearly  seen  that  the  addition  of  the 
optimal  three  runs  has  caused  the  variance  contours  to  have  a  much 
more  circular  appearance,  particularly  in  the  the  case  of  run 
13.  This  demonstrates  the  usefulness  of  the  proposed  measure  as 
an  effective  tool  to  rectify  nonrotatability. 

Example  4.4 

A  company  manufactures  a  liquid  material  for  coating 
automobile  window  glass  which  when  dried  provides  a  barrier  to 
ultraviolet  rays  and  also  reduces  glare  from  the  sun.  The  coating 
liquid  consists  of  water  combined  with  three  active  solid 
ingredients:  a  polymer  (P) ,  a  coupling  agent  (CA),  and  a 

lubricant  (L).  A  film  is  formed  by  passing  I  lie  liquid  through  an 
extruder,  oven  drying  it  and  then  placing  it  on  a  roller. 

An  experiment  is  to  be  conducted  consisting  of  several 
different  combinations  of  the  amounts  of  P,  CA ,  and  1.  with  a  fixed 
amount  of  water.  The  objective  of  the  experiment  is  to  determine 
the  combination  of  P,  CA ,  and  I,  that  is  most  effective  in  terms  of 


reducing  light  penetration.  Light  penetration  is  measured  by 
taking  a  bright  colored  cloth  affixed  to  the  coated  glass  and 
exposing  it  to  light  for  a  fixed  period  of  time.  The  color  of  the 
cloth  is  then  compared  to  that  of  an  unexposed  piece  of  cloth  of 
the  same  color  and  recording  the  degree  of  fading  that  has 
occurred.  A  low  percent  fade  value  is  considered  to  be  desirable. 

In  making  up  the  combinations  of  water,  polymer  (P),  coupling 
agent  (CA),  and  lubricant  (L)  to  produce  the  different  coatings,  a 
central  composite  rotatable  design  in  P,  CA,  and  L  is  set  up. 
Listed  in  Table  5  are  the  amounts  (in  grams)  of  P,  CA,  and  L  to  be 
combined  with  2500  ml  of  water.  The  following  transformation  was 
applied  so  that  the  desing  settings  in  the  factorial  portion  have 
the  familiar  ±1  values: 

P-250  CA-22.5  L-7.5  ,,  ^ 

X1  =  ” 25  ’  *2‘-2.5  ■  *3  ‘-2^-  (4'2) 

The  design  settings  in  terms  of  values  of  Xj ,  X£,  and  x^  are 
displayed  in  Table  6. 

In  this  experiment,  high  ratios  of  the  total  amount  of  the 
active  solid  ingredients  to  water  are  considered  undesirable 
because  of  problems  with  solubility.  For  the  amount  of  water  used 
(2500  ml),  it  was  determined  that  a  total  of  more  than  305  gm  of 
solid  ingredients  would  be  undesirable.  From  Table  5  we  note  that 
design  points  No.  8  and  10  fall  in  that  category.  Their  settings 


had  to  be  reduced  to  conform  to  the  305  upper  bound  constraint. 


The  experimenter  chose  to  reduce  the  levels  of  P  for  these  points 
from  275  and  292.05  to  262  and  275,  respectively.  The  levels  of 
CA  and  L  were  not  altered.  The  total  amounts  of  active  solid 
ingredients  for  these  two  points  are  now  297  and  305  gm 
respectively.  The  design  settings  in  terms  of  values  of  Xj ,  X2 , 
and  for  the  resulting  design  are  given  in  Table  6.  The 
corresponding  percent  rotatability  value  is  81.69,  a  drop  of 
18.31%  from  the  original  central  composite  design.  This  clearly 
demonstrates  the  impact  that  this  change  of  design  settings  has 
had  on  rotatability. 

To  recover  the  loss  of  rotatability,  repair  points  were  added 
to  the  central  composite  design  with  points  No.  8  and  10  modified 
as  was  described  earlier.  The  experimental  region  in  the  space  of 
the  Xj,  X2,  and  X3  variables  is  a  sphere  of  radius  / 3  centered  at 
the  origin.  The  first  optimal  point  selected  by  the  computer 
program  SEARCH  is  shown  in  Table  6  as  point  No.  17.  The 
corresponding  total  of  P,  CA,  and  L  is  256.  / 7  gm,  which  does  not 
exceed  the  305  limit.  This  point  is  therefore  considered 
admissible.  The  percent  rotatability  value  for  the  resulting  17- 
point  design  is  88.79.  This  represents  a  relative  percent 
increase  in  rotatability  of  about  8.7.  The  second  point  selected 
by  SEARCH  increased  the  rotatability  value  to  95.31%.  However,  at 
this  point  x^  =  1.617,  x^  =  .120,  and  x^  =  .119,  which  results  in 
a  total  of  321.023  gn  of  solid  ingredients.  The  point  is 


inadmissible  and  was  subsequently  dropped. 


For  the  next  point, 


the  optimization  was  restricted  to  a  smaller  region,  namely,  a 
sphere  of  radius  /.98  centered  at  the  origin.  The  optimal  point 
(No.  18  in  Table  6)  selected  by  SEARCH  in  this  region  was 
admissible  with  a  corresponding  total  of  304.9  gm  of  solid 
ingredients.  With  the  addition  of  the  latter  point ,  the 
rotatability  value  increased  to  90.83%.  At  this  stage,  the 
sequential  procedure  was  terminated.  It  was  felt  that  additional 
points  inside  the  smaller  region  would  only  produce  marginal 
increases  in  the  percent  rotatability.  To  rapidly  restore 
rotatability,  SEARCH  should  be  allowed  to  operate  near  the 
periphery  of  the  experimental  region  (with  radius  equal  to  / 3  ). 
In  this  example,  however,  this  can  Lead  to  inadmissible  points  as 
was  seen  earlier. 

5.  CONCLUDING  REMARKS 


When  a  design  is  non  rotatable ,  it  may  be  of  interest  to 
assess  the  extent  of  its  departure  from  rotatability.  This  can  be 
particularly  useful  in  situations  where  a  design  is  required  to 
possess  several  desirable  characteristics,  including  rotatability. 
It  is  important  here  not  to  confuse  priorities  when  considering  a 
choice  of  a  response  surface  design.  The  introduction  of  the 
rotatability  measure  does  not  mean  that  rotatability  should  be 
emphasized  at  the  expense  of  other  design  criteria.  On  the 


1 


IV' 

B 


contrary,  the  variance  and  bias  design  criteria,  for  example,  car. 


-22- 


KmJTm  rn.tfm.lf  mjfM.  M.A  A.1  A*.  .A-*  ■  -•*-  lAt  m  *-«-  ■  'SID-h  A  >  JliiAjtAA  »  v  fcVW  kV  JVlfc  X  J If  -VA  A  A 


be  by  far  more  important. 


Only  after  these  criteria  (or  perhaps 


others  that  may  be  of  interest  to  the  experimenter)  have  been 
fully  pursued  might  one  consider  rotatabi 1 i t y.  In  fact, 
rotatability  may  sometimes  be  compromised  in  favor  of  other 
desirable  design  features. 

The  measure  of  rotatability  introduced  in  this  paper  gives 
the  experimenter  greater  freedom  and  flexibility  in  selecting  from 
a  pool  of  efficient  designs  those  that  have  high  percent  rotat¬ 
ability.  These  designs  can  be  made  even  more  rotatable,  if 
necessary,  by  proper  addition  of  design  points  as  was  shown  in 
Section  3.  Certainly,  this  action  will  not  reduce  the  efficiency 
of  the  augmented  designs. 
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APPENDIX  A:  THE  DERIVATION  OF  FORMULA  (2.14) 

The  u>„  vectors  that  appear  in  formula  (2.9)  are  pairwise 
~zm 

orthogonal.  Formula  (2.13)  can  therefore  be  written  as 


d  d 

Q  (D)  =  [  u* '  ( Z '  Z)  -  I  k_  toi  ][u*(Z'Z)~  E  k 
n  rv  L ~  ^  /v  /m  —  2m  —  —  ~ 


2m  ~2iJ 


=  Bu*'  (Z'Z)I!2-  2  I  k.  u* '  ( Z '  Z)w „ 

^  ~  ~  ^  4tn  —  —  —  —2m 

m=2 

d 

+  E  kJ  IIoj„  II2. 

„  2m  ~2m 
m=2 

By  differentiating  QnQ)  with  respect  to  K2m(m=2*  *  *  *  >d)  and 
equating  the  derivatives  to  zero  we  get 
3Q  (D) 

— n— -  =  “  2u* '  ( Z '  Z)o)„  +  2k  II  w-  II2  =  0,  m=2,3,...,d. 

3k„  ~ - 2m  2m  ~2m 

2m 

By  solving  for  <n  we  obtain 


=  u*  '  (  Z '  Z)  u>_  /H  ui.  B2,  m=2 , 3 . d. 

2m  ~  ~  ~2m  —  2m 


The  solution  (A. 3)  satisfies  the  constrain  t  K2m  >  0  ("=2,3 . 

since  u*'(Z'Z)u„  >  0  for  m=2,3,...,d.  This  follows  because 

—  —  —  —2m 

u*'(Z'Z)u>„  is  a  positive  linear  combination  of  the  even  design 
~  ~  ~  ~2m 

moments  of  order  6  =  2ra  (m=2 ,3, . . .  ,d)  for  the  design  J).  It  is 
easy  to  verify  that  the  solution  (A. 3)  represents  an  absolute 
minimum  that  has  the  value 


min[Q  ( D) ]  =  llu*(Z'Z)ll2  -  E„[u*,(Z,Z)tu„  ] 2  / R to  B2  .  (A 

1  n  ~  J  ~  ~  ~  m=2L~  ~  ~  ~2mJ  ~2m 


Table  1.  The  Elements  of  the  Vectors  v(Z'Z),  u(Z'Z),  u*(Z'Z), 

(!>„,  w.  ,  and  the  Diagonal  Elements  of  the  Matrix  £  for 

MJ  ~z  ~  4 

the  32  Factorial  Design  (Example  4.1) 
v(Z'Z)  A  u(Z'Z)  u*  (Z  *  Z)  ojq  w2 


Table  2.  Roquemore's  (1976)  Three-Variable 
Hybrid  Designs  (Example  A. 2)* 


310 

31 1A 

|  31  IB 

X1 

x2 

x3 

X1 

x2 

x3 

X1 

x2 

x3 

0 

0 

1.2906 

0 

0 

2 

0 

0 

61'2 

0 

0 

-.136 

0 

0 

-2 

0 

0 

-61'2 

-1 

-1 

.6386 

_2 1  /  2 

_2 1  /  2 

1 

-.7507 

2.1063 

1 

1 

-1 

.6386 

2 1  /  2 

_2 1/2 

1 

2. 1063 

.7507 

1 

-1 

1 

.6386 

_2 1/2 

2l/2 

1 

.7507 

-2.1063 

1 

1 

1 

.6386 

2 1/2 

2l/2 

1 

-2.1063 

-.7507 

1 

1.1736 

0 

-.9273 

2 

0 

-1 

.7507 

2.1063 

-1 

-1.1736 

0 

-.9273 

1  -2 

i 

0 

-1 

2.1063 

-.7507 

-1 

0  1 

.  1736 

-.9273 

0 

2 

-1 

-.7507 

-2.1063 

-1 

0  -1 

.  1736 

-.9273 

0 

-2 

-1 

-2.1063 

.7507 

-1 

0 

0 

0 

0 

0 

0 

*  The  first  digit  in  a  design  title  is  the  number  of  variables 
the  next  two  digits  are  the  number  of  points.  A  letter 
differentiates  designs  of  the  same  size. 


-2  7- 


Table  4.  Repair  Points  and  Percent  Rotatability  Values 
for  the  Corresponding  Augmented  Designs  for  Example  4.3 


a 


VS* 

o 


m 

l > 


(The  Initial  Design  is  Given  in  Table  3  and  lias  a 
Percent  Rotatability  Value  of  80.65) 


Run  x. 


r-2  Augmented  Design  %  Rotatability  Increase  in 

%  Rotatability 


11  -.1188  -1.8593  D(1)=  D(0)  Plus 


Run  No.  11 


2  -.8295  .0091  D(2)=  D(1)  P 


Run  No.  12 


13  -.1450  -.2764  D(3)=  D(2)  Plus 


Run  No.  13 


89.99 


97.03 


16.38 


S::> 

Run  No. 

P 

CA 

L 

Total  Amount  of 

"u* 

Active  Solids 

TABLE  6.  The  Rotatable  Central  Composite  Design  Unaer  the  Coding 


of  Eqs.  (4.2),  the  Design  Settings  for  the  Modified 
Design  with  the  Additional  Repair  Points,  and  the 
Corresponding  Rotatability  Values  (Example  4.4) 


Rotatable  Central  Modified  Central 
Run  No. _ Composite  Design _ Composite  Design  %  Rotatability 


X1 

x2 

x3 

X1 

x2 

x3 

1 

-1 

-1 

-1 

-1 

-1 

-1 

2 

1 

-1 

-1 

1 

-1 

-1 

3 

-1 

1 

-1 

-1 

1 

-1 

4 

1 

1 

-1 

; 

1 

1 

-1 

5 

-1 

-1 

1 

-1 

-1 

1 

6 

1 

-1 

i 

1 

-1 

1 

7 

-1 

1 

1 

-1 

1 

1 

8 

1 

1 

1 

.48 

1 

1 

9 

-1. 

682 

0 

0 

-1.682 

0 

0 

10 

1. 

682 

0 

0 

1 

0 

0 

11 

0 

-1. 

682 

0 

0 

-1.682 

0 

12 

0 

1. 

682 

0 

0 

1.682 

0 

13 

0 

0 

-1. 

682 

0 

0 

-1.682 

14 

0 

0 

1. 

682 

0 

0 

1.682 

15 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

i 

0 

0 

0 

81.69 

17* 

-.828 

-.506 

-.506 

88.79 

18* 

.966 

.151 

.151 

90.83 

♦Repair  points 
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